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Abstract 

We study a model Dyson-Schwinger equation for the quark propagator closed 
using an Ansatz for the gluon propagator of the form D(q) ~ q 2 /[(q 2 ) 2 + 6 4 ] 
and two Ansatze for the quark-gluon vertex: the minimal Ball-Chiu and the 
modified form suggested by Curtis and Pennington. Using the quark conden- 
sate as an order parameter, we find that there is a critical value of b = b c 
such that the model does not support dynamical chiral symmetry breaking 
for b > b c . We discuss and apply a confinement test which suggests that, for 
all values of b, the quark propagator in the model is not confining. Together 
these results suggest that this Ansatz for the gluon propagator is inadequate 
as a model since it does not yield the expected behaviour of QCD. 

Pacs Numbers: 12.38.Aw, 12.38.Lg, 12.50.Lr, 12.90.+b 
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I. INTRODUCTION 



Dynamical Chiral Symmetry Breaking (DCSB) and confinement are two crucial features 
of quantum chromodynamics (QCD). Indeed, it might be argued that a realistic model of 
QCD should manifest both of these properties since they are responsible for the nature of 
the hadronic spectrum; DCSB ensuring the absence of low mass scalar partners of the pion 
and confinement ensuring the absence of free quarks, for example. 

A natural method for studying both DCSB and confinement in QCD, and models thereof, 
is the complex of Dyson-Schwinger Equations (DSEs). U The equations for the two-point 
functions of gluons and quarks have been used in many such studies. This manifestly 
relativistically covariant approach, recent reviews of which can be found in Refs. has 
provided the foundation for a useful and successful understanding of the phenomena of low 
energy QCD by facilitating the construction of realistic field theoretic models. Q 

Another goal in the DSE studies is to develop this nonperturbative approach to a point 
where it is as firmly founded as lattice QCD and calculationally competitive. Although 
much remains to be done in order to achieve this goal there has been a good deal of progress, 
especially in the study of Abelian gauge theories where direct and meaningful comparisons 
can be made, and agreement obtained, between the results of lattice and DSE studies. || 

In considering DSE studies it is important to note that they are hampered by the fact 
that there is an infinite tower of coupled equations: the equations for the two-point func- 
tions couple them to other two- and three-point functions; the equations for the three-point 
functions couple them to the four-point and higher n-point functions; etc. A commonly 
used resolution of this problem is to truncate the system at a finite number of coupled 
equations by making Ansdtze for the higher n-point functions. For example, one may study 
the DSE for the quark propagator alone by choosing Ansdtze for the gluon propagator and 
quark-gluon vertex thus closing the system. This is the nature of our study. 

Herein we study the fermion DSE obtained with a model gluon propagator (2-point 
function) which vanishes at q 2 = 0: 

1 q 2 
= g 2 [l-II(g 2 )] = (q 2 ) 2 + ' (1) 

where IT(g 2 ) is the gluon vacuum polarisation function and b is a real parameter, in order 
to determine whether it can support DCSB and/or generate a confining quark propagator. 
Two properties we require of the propagator of a confined particle are: 1) the absence of a 
Kallen-Lehmann representation and 2) no singularity on the timelike q 2 axis. The gluon 
propagator obtained from Eq. (|I|) satisfies both of these requirements. 

Such a form for the gluon propagator, even though it may be argued to describe a con- 
fined gluon, is perhaps counterintuitive, since it would appear to provide a weak interaction 
between quarks at small p 2 , which corresponds to large distances. Indeed, some studies of 
the fermion DSE in QCD have employed quite a different Ansatz: one which behaves as 
l/(q 2 ) 2 for q 2 ~ 0. |7]|| This form for the infrared (IR) behaviour of the gluon propagator is 
suggested by a number of studies of the DSE for the gluon propagator in both axial |J and 
covariant gauges [ p!U| , pT| using various approximation and/or truncation procedures; notably, 
they all effectively neglect the 4-gluon vertex. (We note that the results of the axial gauge 
studies |S| may be questioned on the basis that the gluon propagator therein is inconsistent 
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with the known spectral representation in axial gauge. [O) In addition, this form of gluon 



propagator in the infrared is consistent with area law behaviour of the Wilson loop |13 
which has been observed in lattice gauge theory studies of QCD |14| and is often regarded 
as indicating confinement. 



However, the form in Eq. ([I]) is suggested by a number of studies. It has been argued (15 
that in order to completely eliminate Gribov copies, JT6J and hence to fix Landau or Coulomb 
gauge uniquely in lattice studies, one must introduce new ghost fields into QCD in addition 
to those associated with the Fadde'ev-Popov determinant in the continuum. Analysing 
the lattice action thus obtained suggests that the gluon propagator vanishes as (g 2 ) 7 , with 
7 > not determined. Subsequent analysis of a simplified model yields 7 = 1 and, in fact, 
a gluon propagator of the form in Eq. (|l]) with b a finite constant in Landau gauge. Similar 
considerations in Ref. |16[ yield the same result. A propagator of the form in Eq. ([!]) was 



also suggested in Refs. |17H19| as a result of an analysis of an approximate DSE for the gluon 
propagator. A recent lattice QCD simulation [^D] also provides some support for this form. 
We also note that a gluon propagator of the type in Eq. (HD may arise in the field strength 
approach to QCD. jn|] 

It is therefore important to study the phenomenological implications of Eq. ([!]); i.e, to 
determine whether it can support DCSB and confinement in QCD based models of the 
type in Refs. ||. It has been argued |i5||iT| -|i9H that Eq. (|T|) represents confined gluons 
because there are no poles on the timelike real axis in the complex-g 2 plane and it allows the 
interpretation of the gluon as an unstable excitation which fragments into hadrons before 
observation (in a time of the order of 1/6). It is also argued [[OJ that such a gluon propagator 
should lead to a quark propagator with similar structure in the complex plane, and hence 
a similar interpretation, but this result has not been proven. Learning sufficient about the 
quark propagator to make inferences about its analytic structure is therefore an important 
part of our study. 

We study the implications of Eq. ([!]) for the structure of the quark propagator and DCSB 
using the fermion DSE. A similar study is undertaken in Ref. (22fl . Our DSE is closed using 
two Ansatze for the quark-gluon vertex: the minimal Ball-Chiu vertex and that of Curtis 
and Pennington. [^] These vertices are free of kinematic light-cone singularities. We find 
that in either case there are regions of DCSB and unbroken chiral symmetry characterised 
by a two-dimensional phase diagram in (5 2 ,lnr) space, where b appears in Eq. ([I]) above 
and t will be introduced below in connection with an ultraviolet modification of Eq. ([[]). 
The phase transition is second order. We also employ a confinement test which suggests 
that, with either of the vertex Ansatze and for all values of b and r, Eq. (|]) leads to a quark 
propagator that is not confining. Given these results it appears that the gluon propagator 
of Eq. ([D is inadequate as as a model, at least in our DSE framework, since it does not yield 
the behaviour expected in QCD. 

In Sec. II we present a model quark DSE which we solve numerically. We also discuss the 
gluon propagator and quark-gluon vertex Ansatze in some detail. In Sec. Ill we evaluate the 
quark condensate, which we use as an order parameter for DCSB, and determine the charac- 
teristic properties of the phase transition. We also demonstrate that the quark propagator 
is not confining. In Sec. IV we discuss and summarise our results and conclusions. 
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II. MODEL DYSON-SCHWINGER EQUATION 



In Minkowski space, with metric g^ u = diag(l, —1, —1, —1) and in a general covariant 
gauge, the inverse of the dressed quark propagator can be written as 

S-\p) =^-m-S(p) = Z-\p 2 ) (tf-M(p 2 )) = A{p 2 ) 0-B{p 2 ), (2) 

with m the renormalised explicit chiral symmetry breaking mass (if present), S(p) the self- 
energy, M(p 2 ) = B(p 2 )/A(p 2 ) the dynamical quark mass function and A(p 2 ) = Z^ 1 {p 2 ) the 
momentum-dependent renormalisation of the quark wavefunction. 
The unrenormalised DSE for the inverse propagator is 

S-\p) ^-m^-i^f -0^rS(k)T»(k,p)D, u ((p-k) 2 ), (3) 

where Y u is the proper quark-gluon vertex and D^ u (q 2 ) is the dressed gluon propagator. 
Hereafter we set m bare = so that we can concentrate on dynamical symmetry breaking 
effects. (A study of explicit chiral symmetry breaking in such DSE models of QCD can be 
found in Ref. |25| , |2"6|| . ) The renormalised, massless DSE is 

S R \p) = Z s i,- % Z T \g 2 J A ^- 4 rS R (k)TUk,p)D^((p - k)% (4) 

where Zs and Z? are quark-propagator and quark-gluon-vertex renormalisation constants, 
respectively, which depend on the renormalisation scale, fi, and ultraviolet cutoff, A. Ne- 
glecting ghost fields and explicit 3-gluon vertices, as we do in the following (see Eq. (^) 
below), one has Z s = 1 = Z r at one- loop in Landau gauge. Using this result here leads to 
our model renormalised DSE for the quark self energy: 

S(p) =tlg 2 f -0^rS(k)T»(k,p)D, u (( P -k) 2 ), (5) 
where here and hereafter we suppress the label R. 



A. Model Gluon Propagator 



In a general covariant gauge the dressed gluon propagator, which is diagonal in colour 
space, can be written: 



D» v {q 2 



sT- 



q 2 



U(q 2 



q^q v 
q z 



(6) 



where II(g 2 ) is the gluon vacuum polarisation and a is the gauge parameter. It should be 
noted that in covariant gauges the longitudinal piece of this propagator is not modified by 
interactions. This follows from the Slavnov- Taylor indentities in QCD. (See, for example, 
Ref. [|27|, pp. 42-45.) 

The form of the propagator we consider herein has been argued to arise in Landau gauge 
(a = 0) and we shall use that gauge hereafter. (This has the benefit that the gluon propagator 
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is purely transverse. In addition, Landau gauge is a fixed point of the renormalisation group 



and therefore the gauge parameter does not "run". [See, for example, Ref. |27| , pp. 135- 
136.]) 

A commonly used approximation is to write 

where a(—q 2 ) is the running coupling in the gauge theory. In QCD this amounts to neglecting 
ghost contributions to the gluon vacuum polarisation. [ll| (In Abelian theories, of course, 



this is not an approximation but an exact result which follows from the fact that, in the 
absence of ghost fields and gauge-boson self-interactions, the renormalisation constants for 
Il(g 2 ) and g 2 are the same. For Q 2 = —q 2 ^> Aq CD , the one loop contribution is 



dominant and the running coupling is 

, 12 

(We use Nf = 4 herein. Our results are insensitive to changes in this value.) In QCD, 
using Eq. (|7|) in Eq. @ is expected to be accurate at large Q 2 where it leads to the correct 
perturbative leading-log behaviour. |29| 

In considering the small q 2 behaviour we return to Eq. (|l]). We observe that in 
Refs. pl7|-|T9H a solution of the coupled DSEs for the gluon propagator and triple-gluon 



vertex was obtained, using rational polynomial Ansatze, (the coupling to the 4-gluon vertex 
was eliminated) and it was argued that Eq. ([[]) is a good approximation to the solution. 
A significant result in this study is that, in order to support a solution of the type in 
Eq. ([I]), the transverse part of the triple-gluon vertex, Tst(pi,P2,P3), necessarily has kine- 
matic light-cone singularities of the form l/(p 2 ). These singularities, which cannot arise in 
perturbation theory, are introduced in "solving" the Slavnov- Taylor identities for the triple- 



gluon vertex. In this connection we remark that Ref. |23| might be argued to suggest that 
the Slavnov- Taylor identities for both the triple-gluon vertex and quark-gluon vertex can be 
solved without introducing kinematic light-cone singularities. This latter approach would 
allow for nonperturbative corrections to the vertices and a simple connection with the per- 
turbative vertices in the asymptotically free region. In this case the light-cone singularities 



in the triple-gluon vertex of Refs. [p78| , p~9[1 might be viewed as undesirable 



We note too that there has been an attempt to employ the gluon propagator described by 
Eq. ([]]) in a study of quarkonium spectra. |!0| A Blankenbecler-Sugar reduction of a ladder- 
like approximation to the Bethe-Salpeter equation is used and it is argued that in the bound 
state equation that results one can approximate the effect of Eq. (|l|) by a Coulomb poten- 
tial for all r. The interesting structure in this equation arises because the bound fermions 
are described by propagators with complex conjugate poles just as the gluon propagator 
has (see Eq. ([40]) and associated discussion below). The results in this study, however, are 
not competitive with detailed fits using potential models. It may be that the approxima- 



tions/truncations made in Ref. [30] are partly responsible for this. 

In connection with the small q 2 behaviour of the quark-quark interaction the lattice 



results of Ref. |2(| are also interesting. This study, on 16 3 x 40 and 24 3 x 40 lattices at 



(3 = 6.0, yielded a gluon propagator in a lattice Landau gauge which allowed a fit of the 
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form in Eq. (P at small q 2 but which could not rule out a fit using a standard massive 
particle propagator. Other lattice sizes and values of (3 were also studied. The results at 
(3 = 6.3 on a lattice of dimension 24 4 were not inconsistent with these results but in this case 
the small physical size of the lattice was a problem. On a lattice of dimension 16 3 x 24 at 
/? = 5.7 it was found that the gluon propagator was best fit with a standard massive vector 
boson propagator with mass ~ 600 MeV. (We note that the gauge fixing in this study did not 
include the modifications suggested by Ref. [JL5|J.) These studies represent an improvement 



in both technique and lattice sizes over earlier lattice studies of the gluon propagator |31 
but the conclusions are not markedly different. The studies of Ref. [ ffl]| , using (3 = 5.6,6.0 
on a 4 3 x 8 lattice and (3 = 5.8 on a 4 3 x 10 lattice, obtained results that were consistent 
with a free massive boson propagator with mass ~ 600 MeV. Clearly, further lattice studies 
would be of great interest. 

At present, one can only say that Refs. [p~5| — pO|l suggest that in QCD 

1 - (?2)2 (9) 



1 - Tl(q 2 ) (q 2 ) 2 + b A 

is not implausible, at least at small q 2 . For this reason we study Eq. (|J) in the DSE approach 
in order to determine whether it can lead to DCSB and a confining quark propagator. This 
will provide further insight into the validity of this form of gluon propagator. 

Combining Eqs. (|7|) and (^) leads to the ultraviolet-improved model gluon propagator 
that we consider herein; i.e, we study the phenomenological implications of a model gluon 
propagator obtained with 



9 2 1 , (I 2 " 2 



a ^ -1 ) 7^V2 T U » ( 10 ) 



4tt l-n(g 2 ) v ' ^ ' (g 2 ) 2 + 6 4 



where (Q 2 = — q £ 



a(r; Q 2 



dir 



In 



T 



Q 2 



A 2 

QCD 



(11) 



in Eq. @. Here r > 1 is an IR regularisation parameter introduced so that the logarithmic 
singularity is shifted to Q 2 = — tAq CD which ensures that the piece derived from Eq. (g) 



dominates in the spacelike IR region. p2| . |33 



B. Model Quark-gluon Vertex 

In choosing an Ansatz for the vertex we note that in QCD the Slavnov- Taylor identity |54j 



q^(k,p) [1 + b(q 2 )} = [1 - B(q,p)} S'^k) - S'^p) [1 - B(q,p)} , (12) 

where q — (p — k), b(q 2 ) is the ghost self energy and B(p,q) is the ghost-quark scattering 
kernel, constrains the longitudinal part of the vertex. Clearly, the often used "rainbow" 
approximation, in which T^k^p) is simply taken to be 7^, cannot satisfy Eq. ( |12"D in an 
interacting theory. 
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Neglecting ghosts, as we have done so far, this relation takes the form of the Ward- 
Takahashi identity in QED: 



q^(k,p) = S- 1 (k)-S- 1 (p) 



(13) 



The constraints that this relation place on the vertex can now be inferred from the QED 
studies. Taking these into account one is lead to a vertex of the form 



where 



T Bc(k,p) 

T cp(k,p) 
with d(k,p) 



Atf) + A(k 2 ) _ M 

+^{[W-M*)]^-[B<P<)-Bi* 

Y{k 2 -p 2 ) - (k+p) v (f(- f) 



(k 2 +p 2 



2d(k,p) 
(k 2 -p 2 ) 2 + 



[A(k 2 ) - A(p 2 



B 2 (k 2 ) B\p 



+ 



A 2 (k 2 ) A\p 2 



(14) 
(15) 

(16) 
(17) 



In these equations, Tbc is the Ball-Chiu vertex of Ref. [E3[ and Tcp is the additional piece 



suggested and studied by Curtis and Pennington in Refs. [24 



Equation flT4] ) specifies our vertex Ansatz which has the properties that it satisfies the 
Ward-Takahashi Identity, is free of kinematic singularities, reduces to the bare vertex in the 
absence of interactions, transforms correctly under charge conjugation and Lorentz trans- 
formations and preserves multiplicative renormalisability in the quark DSE. (Of these prop- 
erties the minimal Ball-Chiu vertex satisfies all but the last.) 

For the most part in the following we neglected TQ P (k,p); i.e, we used a minimal Ball- 
Chiu Ansatz. As we show below, this has the virtue of simplifying the integral equations. 
We did include this term for a single value of r and a number of values of b in Eq. ( |TUD 
and found, as shown below, that it generates a small quantitative change in some of the 
characteristic quantities calculated in the model but does not alter its qualitative features. 



C. Wick Rotation and Euclidean Space 

The model is now almost completely specified. Up to this point, however, we have not 
considered the possible complications that may arise in employing a Wick rotation to obtain 
the DSE in Euclidean space where a numerical solution is most easily obtained. In all model 
fermion DSEs that have been studied so far the Wick rotation is not allowed in the sense 
that the rotation of the ko contour always encounters at least a pole and very often a branch 
cut. P2"| , PB]| Indeed, in some models it is not possible to rotate the contour at all. |37J This 



point is discussed in Refs. We shall simply complete the definition of our model DSE by 
specifying it in Euclidean space, with metric 5^ = diag(l, 1, 1, 1) and with 7 M hermitian: 

S(p) = lg 2 f -0^rS(k)T»(k,p)D, u ((p-k) 2 ) (18) 
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where 



S 1 (p) = Z7 • p + E(p) = «7 • pA(p 2 ) + -B(p 2 



(19) 



and all of the other elements in this equation are taken to be specified by the expressions 
given above evaluated at Euclidean (spacelike) values of their arguments. 

Using Eqs. (|B| |ID|,|H| , |I5| , |IB| ) we obtain the following pair of coupled integral equations for 
the scalar functions that specify the model quark propagator: 



A(p 2 



B{p 2 



1 + 



16tt f A d A k a(r;(p-kf 



(2tt) 4 (p-k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 



A(k*) A{k2) t A{p2) - 2 [3p-k-h(p,k)) 



2 

2 



-A(k 2 )AA(k 2 ,p 



p 

„2 ( k ■ P) 



P 



- + —n-h(p, k) 



-B(k 2 )/\B(k\p 2 ^ V ' k) 



P 



IQtt 



A Ji 



d 4 k a(r;(p-k) 2 ) 



1 



3 ./ (2tt) 4 (p-k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 
3B(k' J 



x 



L D n^Mk 2 )+A(p 2 ) 



[ v ' 2 
+ [B(k 2 )AA(k 2 ,p 2 ) - A(k 2 )AB(k 2 ,p 2 )] h(p,k) 



(20) 



(21) 



where h(p, k) = 2 [k 2 p 2 - (k ■ pf] /q 2 and AF(k,p) = [F(k 2 ) - F(p 2 )]/[k 2 - p 2 ]. 

Including the additional Curtis-Pennington term in the vertex, Eq. ([16]), these equations 
are modified as follows: 



A(p 2 ) = RHS of (20) 



+ 



16n 



d A k a(r;(p-k) 2 ) A(k 2 )AA(k 2 ,p 2 ) (k 2 - p 2 ) 3(k 2 - p 2 )k ■ p 



3 J (2tt) 4 (p-k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 2d(k,p) 
B(p 2 ) = RHS of © 



pf 



(22) 



(23) 



+ 



16tt 



d A k a(r-(p-k) 2 ) B(k 2 )AA(k 2 } p 2 ) (k 2 - p 2 \ 
(27) 1 (p - k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 2d(k,p) 



■3(k 2 -p 2 



These equations were solved numerically by iteration on a logarithmic grid of 
x = P 2 /Aq CD and y = A; 2 /Aq CD points. In doing this we ensured that our results were 
independent of the seed-solution and grid choice. Our results were independent of the UV 
cutoff, which was A 2 = 5 x 10 8 Aq CD , and this value was also sufficient to ensure that the 
leading- log behaviour of the mass function, Eq. (^) below, had become evident. 



III. ANALYSIS OF ORDER PARAMETERS 



S 



A. Quark Condensate 



We are interested in determining whether the model gluon propagator specified by 
Eqs. (H) and fliTf ) can support DCSB - a crucial feature of QCD. The quark condensate, 
which is gauge-invariant, is an order parameter for DCSB and it is easily related to the 
trace of the quark propagator, which is the focus of our DSE study: 

<5?> " = "J? k ( a£) (a£) dS ' sMs^Bisf j ■ < 24 > 

where d = 12/(33 — 2n/) and \x is the renormalisation point for the condensate, which is 
usually fixed at 1 GeV 2 . This is the parameter that is used to study DCSB in lattice QCD 
and in many of the model studies in the continuum. (It is clear from Eq. (^4j) that an 
equivalent order parameter is B(s = 0) since if this is zero then so is the condensate.) 

As we remarked above, Eqs. fllQl) and ([TT|) in the quark DSE ensure that the leading-log 



behaviour of QCD is retained in the model so that 



2|a _^(H^ cd] )-^ 

This provides another means of extracting the condensate and hence a check on its evalua- 
tion. 

In cases for which our iterative solution procedure for the DSE converged quickly, with 
relative errors of less than 1 x 10~ 6 , the condensate could be obtained easily. However, for 
values of b 2 near a phase transition the convergence could be extremely slow. In those cases, 
the numerical solution was examined at constant intervals through the run (say, every 50th 
cycle), and the condensate evaluated in each case. Aitken extrapolation |38|] was then used 



to find the "infinite-cycles" limit. In several cases the program was subsequently run until 
the solutions had converged to within 1 x 10~ 6 and the extrapolated result always matched 
the actual result to within a few parts in 10~ 6 . 



B. Critical Behavior of the Condensate 

We solved Eqs. (|20|) and (|2~ID for values of lnr in the domain [0.0, 0.7] and b 2 in [0.1, 1.0] 
using the minimal Ball-Chiu vertex and we plot the condensate obtained from our solutions 
in Fig. [I]. This figure shows regions of unbroken and dynamically broken chiral symmetry. 

Our numerical results suggest that the condensate rises continuously from the transition 
boundary and hence that the transition is second order. As a consequence we assumed that 
the order parameter, (qq)^, behaves as 

(qq),(z)^c(l-^y (26) 

for z — > z~ (for z equal to either lnr or b 2 ) and extracted the critical points, z c , and critical 
exponents, /3, using ratio-of-logs methods adapted from Gaunt and Guttmann. |39| We list 



these quantities in Table. | and, in Fig. ^|, plot the critical curve in the (6 2 ,lnr) plane. 
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We also solved Eqs. ( p0[) and ( pT|) with the Curtis-Pennington additions, Eqs. (p2|) and 
(P3f), using lnr = 0.6. The critical curve (in b 2 ) in this case is shown in Fig. |3] along with the 
minimal Ball-Chiu results for the same value of lnr. The effect of the Curtis-Pennington 
addition is to lower the critical value of b 2 but, as we show below, the critical exponent is 
unchanged. This curve illustrates the point that the qualitative features of the model are 
not affected by this modification of the model quark-gluon vertex. 



1. Critical Parameters for the Chiral Phase Transition 
From Table | we find: 

p BC = 0.575 , 073 = 0.024 . (27) 

We note that the critical exponent obtained with In r = is quite different from the others. 
This is a special case since for this value the propagator does not vanish in the infrared: 

^D(q 2 ) = dn^. (28) 

If we neglect this point in our analysis then we find 

p BC = 0.572 , o-p = 0.020 . (29) 



The results in Eqs. ([Z7|) and (^) are in agreement with those of Ref. |33| where it 
is argued that ft = 0.589 ±0.031. That study used b 2 = and found a critical value of 
lnr = 1.69 which complements the results reported herein, as will be seen in Fig. |2|. 

We also calculated the critical exponent using our numerical DSE solutions obtained 
with the Curtis-Pennington addition to the vertex at lnr = 0.6: 

ft cp = 0.579 , ap = 0.015 . (30) 

This suggests that the vertex modification does not alter the critical exponent of the transi- 
tion; a conclusion that is also supported by the observation that the vertex used in Ref. BB 



was not of either of the above forms but was, effectively, a simple modified rainbow approx- 
imation: 

T,(k,p) = A(k 2 ) llM . (31) 



C. Confinement Test 

It is important to determine whether the model gluon propagator specified by Eqs. (^|) 
and ( |T0| ) leads to quark confinement; i.e., the absence of free quarks in the QCD spectrum. 
This form of gluon propagator has been constructed so that the dominant IR behaviour 
ensures that it vanishes at q 2 = 0. In order to determine whether the quark propagator 
obtained as a solution to our DSE can represent a confined particle we follow Ref. and 
adapt a method commonly used in lattice QCD to estimate bound state masses. 
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We write 



2n _ B (P 2 ) 

as{p ] ~ WWY + bW? (32) 



and define 



A * (T ' f) = / ^ liPiT+fS) ^(P 2 ) ■ (33) 

This is the scalar part of the Schwinger function of the model quark propagator. If we now 
define 

A 5 (T) = J d 3 xA s (T,x) (34) 
and, for notational convenience, 

E(T) = -\n[A s (T)] (35) 

then it follows that if there is a stable asymptotic state with the quantum numbers of this 
Schwinger function then 

lim — — — = m ; (36) 

where m > is the mass of this excitation; i.e, this limit yields the dynamically generated 
quark mass. A finite value of m indicates that the quarks are not confined since it ensures 
that the cluster decomposition property is satisfied by this Schwinger function. [|J 

We argue, therefore, that if the limit in Eq. fl36|) exists for a given propagator then the 
associated excitation is not confined: this is our definition of a non-confining propagator. 

In order to illustrate this point we note that the calculation of the "constituent quark 
mass" in the Nambu-Jona-Lasinio model ||T| can be understood in just this fashion: in this 
model m, as define above, is finite and quarks are not confined. 

In contrast, one can consider the model of Ref. PTfl. Applied to this model one finds 



dE(T) t^oo 

~W ~ KT (37) 

where k is a constant and hence the limit in Eq. ( j36|) does not exist. In this case the 
Schwinger function does not satisfy the cluster decomposition property and hence the quarks 
are confined. Alternatively one may say that through self interaction the quark acquires 
an infinite dynamical-mass. This provides another way of understanding the claim that the 
model of Ref. |J7| is confining. 



Another application of this method, which is of direct interest here, is the IR vanishing 
gluon propagator of Eq. ([I]). In this case one has a boson propagator and finds for the 
analogue of Ag(T): 

A(T) oc -i- exp (J?) (cos ( b * )- sin ( b -L) ) . (38) 



11 



We remark that the Schwinger function in this case is not positive definite, which is an 
easily identifiable signal in A(T) that is due to the pair of complex conjugate poles, and this 
violates the axiom of reflection positivity. It follows from this that Eq. (RBI) describes a field 
with a complex mass spectrum and/or residues that are not positive. This is appropriate 
for particles that decay and forms the basis of the argument flT5||IT| -|T9H that a propagator of 



the type in Eq. (|T|) allows coloured states to exist only for a finite time (of the order of 1/6) 
before hadronising; i.e., that the propagator describes confined gluons. 

1. Confinement and Dressed- quark-masses 

In applying this method here it is obvious that numerical evaluation of the Fourier 
transforms required in using Eq. ([!6|) will be hindered by numerical noise as T is increased. 
In order to minimise the effect of this noise on the derivative, we fitted A$(T) to a form 

Cexp(-mT) (39) 

and extracted the derivative from this fit. (Importantly, we found no indication of the struc- 
ture suggested by Eq. (|31f ) in our results.) This was particularly useful with the propagators 
obtained using small values of b which had large dynamical masses (as one would expect 
since the condensate is large in this case) and hence a rapid decline with T. 

We applied the confinement test in the following cases: 1) The propagators obtained 
with lnr = 0.1 and b 2 in the range [0.1, 1.0]; 2) The propagator obtained with lnr = and 
b 2 = 0.35 which yields the largest value of — (qq)^ on the (6 2 ,lnr) domain considered; 3) 
Two propagators obtained with (6 2 ,lnr) = (0.1,0.6) - one using the Ball-Chiu vertex and 
another using the Curtis-Pennington addition. The results obtained by fitting the form 
Eq. (|3~9p to our numerical output are presented in Table. |I[ In Fig. f| we present plots of 
E'(T) for the family of propagators obtained with lnr = 0.1 and this clearly illustrates that 
an unambiguous determination of the dressed-quark-mass is possible in this model. It will 
be observed that the mass decreases with increasing b 2 . This is easily understood in terms of 
the chiral phase transition: as b increases beyond b c there is no DCSB and massless current 
quarks remain massless. Since the behaviour of all the other solutions we obtained was 
qualitatively the same as that described by the results presented in Table. [Tl] and Fig. f| we 
infer that the model considered here does not yield a confining quark propagator. We also 
note that the rainbow approximation studies of the fermion DSE with Eq. ([!]) in Ref. |22" 



which address the question of confinement by a direct continuation to Minkowski momentum 
space, found a quark propagator with a pole at timelike p 2 ; i.e, a non-confining propagator. 

We remark that, within numerical noise, the Curtis-Pennington addition made no dif- 
ference to the dressed-quark-mass value extracted in the cases considered and only slightly 
reduced the normalisation constant C. Clearly, the Curtis-Pennington addition leads only 
to a minor quantitative effect in this part of our study too. 

IV. DISCUSSION AND CONCLUSION 



It has been suggested |1| that the model gluon propagator in Eq. ([[]) would lead to a 



confining quark propagator of the form 
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^- p + Co (40) 

(?7 ■ p + d) (27 • p + c\) 

where c G M and C\ G C are constants; i.e., to a fermion propagator with complex conjugate 
poles just as the gluon propagator has; confinement being realised through the absence of 
poles on the real timelike axis. As Eq. (|38|) shows, such a propagator has a characteristic 
signature in A#(T) which we do not see in Fig. ||. In our model then it is clear that a fermion 
propagator of the type in Eq. ([40]) does not arise. This does not eliminate the possibility 
that it can arise in the approach of Refs. |I7|-0], however, since the rational-polynomial 
Ansdtze employed for the vertex functions therein may lead to a completely different quark- 
gluon vertex to that used here. We simply remark that our results suggest that a fermion 
propagator of the type in Eq. (|4"0| ) cannot arise if the quark-gluon vertex is free of kinematic 
singularities. 

In conclusion, we have studied a model DSE (Dyson-Schwinger equation) for the quark 
propagator using a model gluon propagator that vanishes as q 2 — > 0, Eq. (]!]), and a light- 
cone-regular model quark-gluon vertex, Eqs. (|14"|-iT). This is the first study to analyse 



the phenomenological implications of Eq. ([!]) in the framework of the fermion DSE. Our 
results suggest that this model can only support DCSB for values of In r and b 2 less than 
certain critical values (see Fig. ||) and does not confine quarks. Qualitatively similar results 
are obtained in Ref. As a consequence we believe that this model gluon propagator is 
unlikely to be a useful foundation for a chiral-dynamical model of QCD of the general coupled 
Dyson-Schwinger-Bethe-Salpeter equation type considered in Refs. ||. Indeed, taken in a 
broader context, our results do not support the contention that Eq. ([[]) is the correct form 
of the quark-quark interaction at small q 2 in QCD. 
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FIGURES 



i 



FIG. 1. Criticality plot for —{qq)ji as a function of lnr and b . The condensate, — {qq)fi, is in 
units of MeV, scaled to fj, 2 = lGeV, and b 2 is in units Aq CD ; the gluon regulator r is dimensionless. 

FIG. 2. Critical curve for the phase transition in the (lnr, b 2 ) plane. The asterisk is the result 



extracted from Ref. [33 



FIG. 3. Comparison of the —{qq)jl condensate curves for the minimal Ball-Chiu and Cur- 
tis- Pennington Ansatze for the proper vertex. Both curves have lnr = 0.6. Diamonds, o, connected 
with solid lines are the results from the B-C vertex; plus-signs, +, connected with dashed lines are 
results from the C-P vertex. 

FIG. 4. Dressed-quark-mass curves for the family of propagators with the minimal Ball-Chiu 
vertex and lnr = 0.1. The masses are in units of Aqcd- 
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TABLES 



TABLE I. The critical points and exponents extracted for various values of lnr; the cumulative 



result is /3bc = 0.575, with erg = 0.024; excluding the point with lnr = 0.0, /?bc = 0.572 with 
ap = 0.020. 
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TABLE II. Asymptotic dressed-quark-mass values for the family of propagators with lnr = 0.1, 


the propaj 


;ator which showed maximal DCSB (at lnr = 


= and b 2 = 0.35), and for lnr = 0.6, 


b 2 = 0.1 with both the Ball-Chiu and Curtis- Pennington vertices. 
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Criticality plot, -<qq> 1//3 vs. b 2 , ln(T) 
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Figure 1 



Critical curve for propagator with Ball-Chiu vertex 
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